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Abstract 

The model under consideration is a two-dimensional two-component plasma, 
stable against collapse for the dimensionless coupling constant (3 < 2. The combi- 
nation of a technique of renormalized Mayer expansion with the mapping onto the 
sine-Gordon theory provides the full thermodynamics of the plasma in the whole 
stability range of (3. The explicit forms of the density-fugacity relationship and of 
the specific heat (at constant volume) per particle are presented. 
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1 Introduction 

The model under consideration is the two-dimensional (2D) two-component plasma (TCP), 
i.e., a neutral system of pointlike particles of positive and negative charge, interacting 
through the 2D (logarithmic) Coulomb interaction. Classical equilibrium statistical me- 
chanics is applied. The system is stable against collapse for the dimensionless coupling 
constant (3 (sometimes denoted as F) lower than 2. The absence of some short-range 
interaction, e.g. the hard-core one, prevents to extend the treatment beyond (3 = 2 and 
to reach the Kosterlitz-Thouless phase transition which occurs around /3 = 4. 

In arbitrary dimension, the long-range tail of the Coulomb potential causes screening, 
and thus gives rise to exact constraints for the charge-charge structure function (see review 
|I|]). In the bulk regime, the zeroth- and second-moment conditions for the structure 
function are known exactly 0, 0. The specialization to the 2D TCP brings additional 
sum rules. The knowledge of the exact equation of state fixes via the compressibility 
equation the zeroth moment of the number density correlation function p. In a recent 
paper 0, the second- moment formula for that density-density correlation function has 
been derived by using analogies with critical systems. This result was rederived directly 
by using a renormalized Mayer expansion [0, p in Ref. 0. 

At arbitrary temperature, the 2D TCP is supposed to be in the critical state in terms 
of the induced electrical-field correlations (but not the particle correlations) |Tl|. The 



free energy therefore exhibits a universal finite-size correction predicted by the Conformal 
Invariance theory. 

The 2D TCP is mappable onto the Quantum Field Theory (QFT) models, namely 
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the 2D Thirring model (see e.g. [|T2|) and the 2D Euchdean, or (l+l)-dimensional quan- 



tum, sine-Gordon model |13], |14|, [15]. The integrability of the sine-Gordon theory was 



established and analysed in Refs. [|I6[, |T^, [|T^, |T9[, |2y] and in many others. 

This paper is devoted to the derivation of the thermodynamics of the 2D TCP in the 
whole stability interval < (3 < 2. The aim is reached via a study of the relationship 
between the fugacity z and the total number density of particles n. It is known from a 
scaling argument that the combination n^^~/^/'^^ / z is a function of P only. Up to now, the 
exact form of the n — z relationship was known only in the Debye-Hiickel (3^0 limit. 



n/z — > 2, and at the exactly solvable collapse border j3 = 2 [^, |2^, when y/n/ z oo. 
The little interest in the topic in the past is surprising in view of the fact that in the 
QF theories related to the 2D TCP one-point expectations of local fields are the objects 
of primary importance. In the Conformal Perturbation treatment of a QFT, all multiple 
correlation functions are reducible by the Operator-Product-Expansion |2^ to the one- 
point vacuum expectations which contain all the non-perturbative information about the 
QFT. 

We first investigate the n — z relationship by using the renormalized Mayer expansion 
for many- component fiuids and find out that rS^'^/^"* /z ~ P^^^^x an analytic function 
of solely P, /(/?), whose Taylor expansion around /? = can be constructed systematically 
by evaluating certain convergent integrals of Bessel functions. The mapping onto the 2D 
sine- Gordon theory, whose normalization of the relevant cos- field is consistent with the 
formalism of the renormalized Mayer expansion, then provides the explicit form of the 
function /(/?), checked on a few lower orders of its /^-expansion. The exact density-fugacity 
relationship provides the full thermodynamics of the 2D TCP. 



The paper is organized as follows. Section 2 recapitulates briefly the ordinary Mayer 
formalism for many-component fluids. The procedure of the renormalization of the Mayer 
expansion is described in Section 3. The application to the 2D TCP is presented in Section 
4. The information gained from the mapping onto the 2D sine-Gordon QFT is the subject 
of Section 5. The exphcit form of the speciflc heat (at constant volume) per particle is 
given together with a discussion of the results in Section 6. 



2 Ordinary Mayer expansion for many-component 
fluids 

We flrst consider a general classical multi-component fluid in thermodynamic equilibrium 
at temperature T. The last system is composed of distinct species of particles {a}. The 
particles may be exposed to an external potential u{r,a) and interact with each other 
through the pair potential v{ri,ai\rj, aj) which depends on the mutual distance \ri — r,| 
of particles i,j as well as on their types ai,aj. Vector position Tj of a particle in d- 
dimensional space will be sometimes represented simply by i. In the grand canonical 
formalism, denoting by {/ia} the chemical potentials of species {a}, the grand-canonical 
partition function is deflned by 

N 



1 r ^ 

E{T, V, {/.J) = E n TTT / n ^(^' '^Ol exp 



4 Yl '"ihcri\j,crj) 



(1) 



where the flrst sum runs over all possible species numbers, N — J2a is the total particle 
number, (5 — 1/kT the inverse temperature and z{r, a) — exp\Pii^ — (3u{r, a)] denotes the 
fugacity of cr-particles. In the direct format, i.e. with fugacities {z{r,a)} as controlling 



variables, H is the generator for the particle densities {n{r, a)} in the sense that 

6lnE[z] 



(2) 



6z{r,a) 

The transition to the inverse format, i.e. with densities {n{r,a)} as controlling vari- 
ables, is based on the Legendre transformation 

- (3F[n\ =\nE- f drJ2n{r,a)\nz{r,a) (3) 

which defines the Helmholtz free energy F as the explicit density functional. The sub- 
traction of the one-particle part provides the dimensionless (for notational convenience, 
minus) excess free energy 

A[n] = -(3F[n] + / dr^ [n{r, a) \nn{r, a) - n{r, a)] (4) 



It is easy to show that A[n] is the generator in the following sense 



In 



n{r, a) 



z{r,a 



6A\n] 



5n{r, a) 

Its Mayer diagrammatic representation in density reads 



(5) 



A[n] = I all connected diagrams which consist of > 2 field n{i,ai) — circles 

and f{i, ai\j, aj) — bonds, and are free of connecting circlesj (6) 

(the removal of a connecting circle disconnects the diagram). Here, 

fih (^j) = exp [-(3v{i, ai\j, aj)] - 1 (7) 

is called the Mayer function and, besides the integration over spatial coordinate of a field 
(black) circle, the summation over all cr-states at this vertex is assumed as well. 



3 Renormalized Mayer expansion 

The renormalized Mayer representation of A[n] results from Eq.(^) in two steps 
i) the expansion of each Mayer function in the inverse temperature, 



/(I, a,\2, ^2) = -f3v{l, (Ti|2, (T2) + - [-f3v{l, (Ti|2, a^)]' + ... (8) 



or, graphically. 



/ -Pv ^ ^ ^ ^ 

o o = o o + a JO + . . . (8') 



l,ai 2,0-2 1,0-1 2,0-2 1,0-r ~ ' '2,0-2 



where the factor l/(number of interaction lines)! is automatically assumed; 
ii) the consequent series elimination of two-coordinated field circles between every couple 
of three- or more-coordinated field circles (hereinafter, by coordination of a circle we mean 
its bond-coordination, i.e. the number of bonds meeting at this circle). The renormalized 
i^-bonds are given by 
K 

o------^----— o = 0- — — -o + o- • -0-I-... (9) 

1,0-1 2,0-2 1,0"! 2,(72 l,Cri 2,0-2 

or, algebraically, 

K{l,ai\2,a2) = [-Pv{l,ai\2,a2)]+J2 [ d3 [-Pv{l, ai\3, a^)] ^(3,^3) K{3, a3\2, 

(9') 

The procedure of bond-renormalization transforms the ordinary Mayer representation 
(|) of A into 

00 

A[n]= ..---^ + Do[n]+ ^ D,[n], (10a) 

s=l 



where 



Do 



> + 



A 

/ \ 



t- - - -t 

X + ; ; + 

■ ^ 4- - - -* 



(10b) 



oo 1 N 

= ElT77 E / n[d^n(z,aO][-Ml,ai|2,a2)] 

N=2 0-1...0-JV ''^ i=l 

[-Pv{2, a2\3, as)] . . . [-Pv{N, a;v|l, ^i)] 

is the sum of all unrenormalized ring diagrams (which cannot undertake the renormaliza- 
tion procedure because of the absence of three- or more-coordinated field points) and 



^2^3 = \ all connected diagrams which consist of > 2 field 

S = l [ 

n(z, (Tj)— circles of (bond) coordination > 3 and multiple 
K{i, ai\j, aj)— bonds, and are free of connecting circles > 



;ioc) 



represents the set of all remaining completely renormalized graphs. By multiple K- 
bonds one means the possibility of an arbitrary number of if-bonds between a couple 
of field circles, with the obvious topological factor l/(number of bonds)!. The order of 
s-enumeration is irrelevant, let us say 



e 






etc. 



In accordance with relation (|^), ln[r;,(l, o"i)/z(l, ai)] is expressible in the renormalized 
format as follows 



In 



o- - 

1,0-1 



+ (io(l,o-i) + ds{l,(Ti, 



(12a) 



where (io(l; = 5DQ/5n{l^ ai) can be readily obtained in the form 

(12b) 



rfo(l, (Ti) = i lim [K{1, ai|2, as) + ai|2, da) 



and 



4(l,ai)= ff^ (12c) 
()n(l, (Ti) 

(s = 1, 2, . . .) denotes the whole family of (1, (Ji)-rooted diagrams generated from Dg. To 
get a family, one has to consider the functional dependence of the dressed if-bonds (H) 
on the species densities as well. Since relation (9') yields 

= a,|3, a3)ir(3, a,) (13) 

dn(3,(j3) 

the root (white) circle is generated, besides the field-circle positions, also on i^-bonds, 
causing their correct K — K division. For example, in the case of generator Di drawn in 
(|ll|), one gets 

l,cri 

(il(l,C'"x) = C>^^^'•^/^^v/•^^^ -|- (14) 

4 The 2D TCP 

Let us now concentrate on the neutral 2D TCP of positive (cr = +) and negative {a = — ) 
pointlike unit (g = 1) charges, with the Coulomb interaction energy given by 

- I3v{i,ai\j,aj) = a^aj [-f3v{i, j)] (15a) 
v{t,j) = -ln{\t -j\/L) (15b) 

where the constant L is for simplicity set to unity and j3 now plays the role of the 
dimensionless coupling constant. Considering the regime with homogeneous densities 
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n{r,a) = n^j, the requirement of the charge neutrahty imphes n_|_ = n_ = n/2 (n is the 
total number density of particles). In the fugacity language, the chemical potentials of ± 
charged particles must equal to one another, 2;(r, a) = z{r)\ in the infinite-volume limit, 
z{r) acquires its bulk value z. 

The renormalized bonds, defined by (9'), now take the form 

K{i,ai\i,aj) =aiajK{i,j) (16a) 

where K{i,j) satisfy the relation 

K{1, 2) = [-pv{l, 2)] + / d3 [-pv{l, 3)] n K{3, 2) (16b) 

In the V ^ oo limit, characterized by translationally invariant K{i,j) = K{\i — j\), the 
Fourier transformation of (|16b| ) results in 

K{k) = [-(3v{k)] + 27m [-(3v{k)] K{k) (17) 

Since the Fourier component v{k) = l/fc^, one arrives at 

r d'^k 1 
K(r) = —8 / — — exp(iA;-r) 



= -/3Ko{r^2n(3n) (18) 

where Kq is the modified Bessel function of second kind. The special scaling form of K 
has a fundamental impact on the n- and /^-classification of renormalized diagrams. 
Let us first study the renormalized representation (|l^) of the generator A{n): 



(i) The first term on the rhs of ( |10a| ) 



^ ^ / dld2n(l,ai)f7i(72[-A;(l,2)]n(2,a2) (19) 



2! 

(Jl,(72 = ±l 



is fixed to zero by the charge neutrahty. 

(ii) The second term Dq ( 10b ) is expressible as follows 



oo » N 



Doin) = E ^ / n M^(l,2)] [-f3v{2,3)] 

N=2 •'^ i=l 



-PviN, 1)] 



(20) 



or, equivalent ly, 



Doin) = - dn'Y: n'^""-'^ / [] 2)] [-/5t;(2, 3)] . . . [-^v{N, 1)] (20') 

/ JO jY^2 "'^ i=l 



In the V ^ oo limit, with respect to ( |16b|) , the sum over N on the rhs of (20') is nothing 
but V X limr^o[-^(^) + /5'i^('")] evaluated at n = n' . Using the small-x expansion of Kq{x) 



oo 



^o(x) = - In (-)/o(x) + E ^577^^(^ + 1) 



i=0 



(21) 



where 



h{x) = E 



oo 



i=0 



22.ui)2 ^(x) = -lnr(x) 



is the psi function [ip{l) = — C, C the Euler's constant], one finally gets 



Doin) P 6n 
= -(nlnn-n) + — 



(22) 



(iii) Let the given completely renormalized diagram Dg (s = 1,2, . . .), belonging to the 
sum ( |10c| ), be composed of A^^ skeleton vertices i = 1, . . . , A^^ of coordination z/j > 3 and 



Ls bonds a = 1, . . . , Lg. Dg can be formally expressed as 



Ds[n] =ts E / n n{i,ai)] J| K(ai, |a2, o^aa) 

(T1...CTIV. =±1 ^ i=l a=l 



(23) 



where tg is the topological factor and ai, 02 € {1, . . . , A^,,}, ai < 02, denotes the ordered 
pair of vertices joint by the a-bond. Since, according to (|16a|), i^'(Q;i, cTq, J 0:2, ctq-J = 
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(Taj(Ja2-^«, Ka = K{ai,a2), Eq. (p3D can be rewritten as follows 

Ds[n] =t, j2 y n ^0 <i n ^» (24) 

(Tl...(TjVs=il ^1 = 1 Q! = l 

The set of coordination numbers {ui} is constrained by J2f=i = 2-^^ as every bond is 
shared by just two vertices. For n{i, cij) = n/2, one has 

D.{n)=tj-] y: (25) 

cri...o-jVg=±l j=l ^ i=l a=l 

i.e., Ds{n) 7^ if and only if the coordinations of all vertices {ui} in Ds are even numbers 
(> 4). Let us suppose that this condition is fulfilled. In the K — > 00 limit, due to 
the invariance of the integrated product Yl^Li Ka with respect to a uniform shift in all 
integration variables {i}, one of these variables can be chosen as a reference put at the 
origin 0, with the simultaneous multiplication by volume V", 

„ Ns Ls 

Ds{n) = tsU^'^V / n n (26) 

•' 1=1 a=l 

j = 1, . . . , Ng. The scaling form of Ka = —f3Ko{\a\\^2¥(3n) permits us to perform the n- 
and /^-classification of the integral in (^). Every dressed bond brings the factor —jS 
and enforces the substitution r' = ry/27rj3n which manifests itself as the factor l/(27r/3n) 
for each field-circle integration ~ / rdr. Since there are (A^^ — 1) independent field-circle 
integrations in p6|) we conclude that 

£M = n/?^-^^+i4 (27a) 

where ds is an intensive quantity, the number 

Ds(n = l,p = 1) , , 

ds = ^ ^ (27b) 
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The first nonzero diagram from tlie sketcli (|TT]) is D2- It contributes to tlie (5^ order, witli 



do 



d^r 



2!4! J 2-K 



2!4! 8 



(2^ 



where ( is the Riemann's zeta function (see Appendix). In the next order only diagram 



Dft survives, and 



1 



3!(2!)3J 2n 2tt 

- C(3) 
3!(2!)3 16 ^ 



(see Appendix), etc. 

The above paragraphs (i)-(iii) are summarized by formula 



(29) 



A(n) P Bn 
-Y~ = j{nlnn-n) + — 



+ n ^ 4/3 



Ls-Ns+l 



(30) 



s=l 



Here, {dg} are the numbers yielded by the topology of diagram Dg, nonzero only if the 
coordinations of all vertices are even numbers > 4. 

In order to evaluate \n{ncr/zcr) using (|^), we first recall the well-known formula 



dA{n) 
dn 



Ed 



5A[n] dn(r,a) 
6n{r, a) dn 



(31) 



valid for an arbitrary functional A[n] with n{r, a) substituted by some function of n. In 



the considered case n{r, +) = n(r, —) = n/2, relation (pTD takes the form 



dA{n) V / SA[n] 



dn 



2 \5n{r,+) 



+ 



6A\n] 



(32) 



n{r,a)=n/2 / 



n{r,a)=n/2 ^^{r, —) 

where the reference r-independence of the functional derivatives was assumed. For dia- 
grams in A[n], the direct link between the -|- and — states of the root point is realized 
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through the state {ai — transformation at field vertices. The diagrams are invariant 
with respect to this transformation, so that it holds 



SA\n] 



6n{r, +) 



SA\n] 



6n(r, —) 



n{r,a-)=n/2 



dA{n)/V 
dn 



(33) 



n{r,u)=n/2 

The consequent equalities \n{n^/z) = ln{n^/z) = d[A{n)/V]/dn, with A{n)/V given by 
(PPP, lead after simple algebra to the final result 



n 



1-/3/4 



2/3^/^ exp 
2/3'^/" exp 



(34a) 



2C + In ( - 



TT 



1 + ^^(3) (f)' + C(3) (f\\o{P')p4h) 



The exact solution at the collapse border /? = 2 [^] confirms the expected divergence 
\fnjz oo, and therefore the radius of convergence of the series Y^^=\(^^''~^''^^ds corre- 
sponds to < /3 < 2. In the next section, from the equivalence of the 2D TCP with the 
sine-Gordon model, the explicit form of the rhs of ( p4a| ) is found. 

The internal consistency of the formalism presented in this part is confirmed by the 
correct reproduction of the equation of state for the 2D TCP 0], 



(35) 



where P is the pressure. Writing down InH = fiPV in Eq.(|]) and combining with Eq.(|D, 
one obtains 

A f n\ 

/3P = — + n-nln — (36) 
V \2zJ 



The consequent substitution of ( pOD and (|34a|) immediately leads to the exact Eq. 



13 



5 Mapping onto the sine-Gordon model 

It is well known that the 2D TCP is equivalent to the 2D Euclidean (classical) sin-Gordon 
model |13[-|15[. In particular, the grand partition function (|I|) can be turned into 



/ V(p exp (/ dV L) 



JV(j) exp -/dV UV(j)y 



(37a) 



where 



■2 my 



+ 22; cos i\l 271 (3 (p 



(37b) 



The rescaling of the real field 

L 



/Stt transforms L into the form 



IGtt 
4 



— iy(t)f + 2z cos (f3(j) 



(38a) 
(38b) 



which is more convenient for our purposes. If we identify one of the two spatial coordinates 
with an imaginary time, say y = it, L (P8a|) can be viewed as Lagrangian in the (1 + 1)- 
dimensional Minkowski space. Defining the canonical momentum conjugate to the field 



0, 



U{x) = ^ = -0 



(39) 



the Hamiltonian density of the corresponding quantum (1 + l)-dimensional sine-Gordon 



model reads 



H 



1 



+ (V0)' 



— 2z cos 



(40) 



The thermodynamic properties of the 2D TCP at finite temperature may be therefore 
obtained from the ground-state properties of the corresponding quantum sine-Gordon 
model, and vice versa. 
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The 2D Euclidean sine-Gordon theory (|38D has been intensively studied |T6[|-||2CI||. 



It can be regarded as the Gaussian Conformal FT perturbed by the relevant operator 
2cos(/50). As concerns the renormalization of parameters (3 and only z renormalizes 
(multiplicatively) . Under the proposed normalization of the field cos(/50) in terms of the 
short- distance limit of the two-point correlation function 

{cos{B(t>){x)cosCPct>){y))^]^\x-y\-^~^" as |a;-y|^0 (41) 

the field cos(/90) has the dimension [length] so that the dimension of parameter z 
is [length]^'^^"^, in full agreement with ( p4a| ) z ~ n^"^/^ (the dimension of the density n 
is [length]"^ and (3"^ = /3/4). The discrete symmetry of the theory (|38D (p (p + 27m/ (3 
(n integer) is broken in the region < /3 < 1; one has to consider one of infinitely 
many ground states {|0„)} characterized by (</>)„, = 27Tn/f3, say that with n = 0. The 
underlying sine-Gordon model is integrable along the standard lines of the Bethe ansatz 
technique |T^, |T^. The spectrum of particles involves solitons S, antisolitons S and 



soliton-antisoliton bound states {Bfj = 1, 2, . . . < 1/^}, 



with masses 



e = ^ = ^ (42) 



mj = 2Msin( — ) (43) 



where M is the soliton mass. The dimensionless specific grand potential cj, defined by 

^ = -ilnS (44) 
[and identical with the ground state energy of the quantum (1 -|- l)-dimensional sine- 
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Gordon model] was found in Ref. (Eq. 2.67): 



8sin(7rO 



(45) 



Mass mi, associated with the hghtest Si-bound state, is given by relation (|43|), so that 

(45') 



M2 (-rC 

oj = tan — 

4 I 2 



Within the framework of the normalization (^Tj), the parameter z is related to the soliton 
mass M as follows (see [0, Eqs. 2.12 and 4.1) 



r 






vrr 





M- 



2r(|) 



2-2(3^ 



(46) 



where F stands for the gamma function. Notice a different parameter notation in Refs. 

H, [HI- 



Returning to the /3-parameter [see formulae (|38b| ) , (|42|) ] and applying the transforma- 



tion relation between the Gamma functions 25 



T(x) = -r(l + x) 

X 



the elimination of M from Eq.(|^) implies the explicit z-dependence of u, 

TT/j 



(47) 



1 - - I Z 1-/3/4 



r(i + f) 



1-/3/4 



tan 



2(4-/3) 



1 + 



/3 



2(4-/3) 



7r/3 



1 

2(4-/3) 



r2 



2 ^ 2(4-/3) 



(48) 



Inserting then (— tu) into the generating relation for the particle density 

d{-co) 



n = z- 



dz 



(49) 



one finally obtains 

^1-/3/4 



2 ^ 



/3/4 



2; 



r(i-f) 
r(i + l) 



tan 



7r/3 



2(4-/3) 



1 + 



/3 



2(4-/3) 



1-/3/4 



7r/3 



2(4-/3) 



ir2 



1 _L /3 

2 2(4-/3) 



(50) 
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The rhs of ( pOD predicts the stabihty against collapse for /? < 2 and the divergency of n at 
/3 — >• 2, as it should be. Its expansion around /? = up to the term, accomplished e.g. 
by using the symbolic computer language Mathematica, turns out to be identical to our 
previous result ( |34b| ). This is the strong evidence that formula (|50D reflects adequately the 
diagrammatic series P^''~'^''~^^ds in ( p4a| ) , and represents the exact density-fugacity 
relationship for the 2D TCP. 

We would like to stress that there were attempts in the past to get the thermodynamics 
of the 2D TCP from the mapping onto the sine-Gordon theory, but they failed because 



of an inadequate transfer of the information. We mention one earlier work p6|, the 



result of which was proved to be wrong in report The crucial point here is that 

the dependence of the parameter z on the soliton mass M, relation (^), was determined 
within the field normalization (^Tj), the last being consistent with the formalism of the 
renormalized Mayer expansion (Sections 3 and 4). 

6 Thermodynamics of the 2D TCP 

The simplest way to get the thermodynamics of the 2D TCP is to pass from the grand- 
canonical to the canonical ensemble via the Legendre transformation 

F{T,N)=n + ^N (51) 

where 



n = keT uo{(3,n)V (52a) 
-uj{l3,n) = (l-f)n (52b) 



4 
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and 



(53) 



The dimensionless specific free energy /, defined by / = F/{NkBT), is then written as 



r(i + f) 



R\ I tan 



2(4-/3) 



1 + 



2(4-/3) 



7r/3 



2(4-/3) ^ 



r2 



+ 



(54) 

.2 ' 2(4-/3) J J 

According to the elementary thermodynamics, the internal energy per particle, e = 
{E)/N, and the specific heat (at constant volume) per particle, cy — Cy/N, are given by 



d 



e = 



Cy 



For the specific heat, one gets explicitly 

4 



dl3 



(55a) 
(55b) 



A;b 4 4-/3 ' 16 



(4 - py 



(1) 



/5 



47r2/32 cos(^ 



(4-/3)%in^ (^) 



(56) 



^4-/5; " V8-2/5, 
The series representation of the psi function and of its derivatives is presented in Ap- 
pendix, formula (A7). It is seen that the density-independent cy jks has a regular high- 
temperature expansion in powers of /3: 



(57) 



The expansion oi cy jkn around the collapse /3 = 2 point results in the Laurent series 



(58) 
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The leading term coincides exactly with the conjecture of Hauge and Hemmer |^ [see 
their Eq.(34)] based on an independent-pair approximation. 

In conclusion, it was shown that the renormalized reformulation of the Mayer expan- 
sion in statistical mechanics of the 2D TCP produces a direct link to techniques in the 
related QFT. The obtained results might represent a step towards the complete integra- 
bility of the 2D TCP on higher correlation levels. A further exploration of the intimate 
relationship to recent studies of the sine-Gordon model is needed. 
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Appendix 

We aim at proving the relations: 



I 



2Tr 



(Ala) 
(Alb) 



The Fourier component of Kl{r), denoted as G{U), can be expressed in 2D as follows 



In 


2 ~ 








2 



(A2) 



where Jq is the ordinary Bessel function. Simple algebra yields 

j '^Kl{r)Kl{r) = f dk kG\k) 



dk kG^{k) 



(A3a) 
(A3b) 



27r 27r 

The primitive functions of the functions kG^{k) and kG^{k) are available explicitly e.g. 
with the aid of Mathematica. They converge to at A; — > oo, and their values in the A; — > 
limit imply 



HdkkG^k) = ^-2{PolyLog[3,i] + PolyLog[3,-i]} 
Jo I 

POO 3 

/ dkkG^k) = -C(3) + 3{PolyLog[3,i] + PolyLog[3,-i]} (A4b) 
Jo 4 



(A4a) 



where 



PolyLog [n,a;] = ^ — 
j=i J 



(A5) 
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is the nth polylogarithm function of x. It is straigthforward to show that 



PolyLog[3, i] + PolyLog[3, -i] = — 



The n-th derivative of the psi function is given by 



OO 1 

c(3) - E 



J=i [3 - 



OO -1 



Deriving twice the functional relation |25 



V^(2x) 



ipi^x) +'ip [x + - 



In 2 



at a; = 1/2, we have ^^^^1/2) = 7^/'(2)(i). Thus, 



1 



i,.«(i)._I,«(i).rc(3) 



where we have used the definition of the Riemann's zeta function C(3) 
insertion of the consequent equality 



(A6) 



(A7) 



(A^ 



(A9) 

-^(2)(l)/2. The 



PolyLog[3, i] + PolyLog[3, -i] = -- C(3) 

lo 



(AlO) 



into (A4), and the consideration of Eqs.(A3) proves (Al). 
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